Spatially anisotropic kagome antiferromagnet with Dzyaloshinskii-Moriya interaction 
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We theoretically study the spatially anisotropic spin- 1/2 kagome antiferromagnet with 
Dzyaloshinskii-Moriya (DM) interaction using a renormalization group analysis in the quasi-one- 
dimensional limit. We identify the various temperature and energy scales for ordering in the system. 
For very weak DM interaction, we find a low-temperature spiral phase with the plane of the spiral 
selected by the DM interaction. This phase is similar to a previously identified phase in the absence 
of the DM interaction. However, above a critical DM interaction strength we find a transition to 
a phase with coexisting antiferromagnetic and dimer order, reminiscent of one-dimensional antifer- 
romagnetic systems with a uniform DM interaction. Our results help shed light on the fate of two 
dimensional systems with both strong interactions and significant spin-orbit coupling. 
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I. INTRODUCTION 

In recent years, spin-orbit coupling has become a cen- 
tral theme in condensed matter physics j 1 ^ Much atten- 
tion has been focused on time-reversal invariant topolog- 
ical insulators (TI) which are driven by strong spin-orbit 
coupling in the weak interaction limit.— ~— Of central in- 
terest is the nature of the phases realized when both spin- 
orbit coupling and interactions are significant £ A number 
of studies have investigated the fate of the TI phase in 
two and three dimensions as the strength of the interac- 
tions is increased^— which reveals some possible out- 
comes of the competition between interactions and spin- 
orbit coupling. Under some conditions, quantum spin 
liquid o 10 i 12 i 14 and other exotic states are expected,— 
which suggests interesting connections between TI and 
spin systems J£ 

In this work, we study a highly frustrated spin model- 
the spin- 1/2 kagome antiferromagnet with spatially 
anisotropic exchange constants and spin-orbit coupling. 
We investigate the physics of this system as the spin-orbit 
coupling, realized in the form of Dzyaloshinskii-Moriya 
(DM) interaction ; 1 ^ 18 is increased. The ground state of 
the spatially isotropic spin- 1/2 Heisenberg antiferromag- 
net in the absence of DM interaction has been studied for 
some time^~— Its large classical ground state degener- 
acy was thought to make it an excellent candidate for 
a quantum spin liquid, and the most recent density ma- 
trix renormalization group results seem to indicate this is 
indeed the case^ Exactly solvable spin-3/2 spin liquids 
have also been studied on the isotropic kagome latticed 

However, some kagome antiferromagnets, such as vol- 
borthite, Cu 3 V 2 07(OH) 2 -2H 2 0, are spatially anisotropic 
which allows a quasi-one dimensional approach to be used 
in their studyi^ 3 . - — This approach is especially powerful 
because one-dimensional methods are adept at describ- 
ing strongly correlated physics i 38 ' 39 In this work, we ex- 
tend earlier studies of the DM interaction on the spatially 
isotropic kagome lattice^— to the spatially anisotropic 
case by using a quasi-one dimensional approach. Our 
main result is that we find a transition from a spiral or- 



dered state with weak Neel order perpendicular to the 
spiral 3 ^ to a spiral state with weak in-plane Neel order 
coexisting with dimer order as the strength of the DM 
interaction is tuned up. (See Fig[T]and Figl^l) We argue 
this transition may be in an experimentally accessible pa- 
rameter range and could possibly be tuned by pressure. 

Before we give a detailed account of our theoretical 
results, we first begin with a brief survey of the growing 
number experimental spin-1/2 kagome systems. 



A. Brief survey of related experimental systems 

A major motivation for the experimental and theoret- 
ical study of spin-1/2 kagome antiferromagnets is the 
possibility of spin liquid and other highly fluctuating 
phasesi 31 i 53 The kagome lattice leads to geometrical frus- 
tration which works against an ordered state, as does the 
small size of the spin which enhances quantum fluctua- 
tions. However, there is theoretical evidence that the DM 
interaction, allowed by symmetry on the kagome lattice, 
tends to drive a spin liquid towards a magnetically or- 
dered statei 44 ' 46 In practice, then, it is important to have 
an estimate of the size of the DM interaction in a given 
system. Experimental realities such as defects and lat- 
tice anisotropy may also influence the low-temperature 
phases. It is thus important to include these aspects in 
theoretical studies as well. In this paper, we focus on 
the role of lattice anisotropy in spin-1/2 kagome antifer- 
romagnets with a DM interaction. 

Among the experimentally studied spin-1/2 kagome 
materials, an absence of long-range magnetic or- 
der down to small temperatures (compared to the 
characteristic exchange energy) is common. Im- 
portant examples of spin-1/2 kagome antiferromag- 
nets include herbertsmithite, ZnCu 3 (OH) 6 Cl 2 ,— 
vesignieite BaCu 3 V 2 8 (OH)2 f S^ 3 - volborthite 
Cu 3 V 2 7 (OH) 2 • 2H 2 0^^ MgCu 3 (OH) 6 Cl 2 (a cousin 
of herbertsmithite)^ and Rb 2 Cu 3 SnFi 2 ^ 9 . 

Herbertsmithite has isotropic exchange interactions 
and has received perhaps the largest amount of theo- 
retical and experimental attention among the materials 
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listed above. No magnetic order appears in this system 
down to temperatures of 40 mK (nearest-neighbor ex- 
change J ~ 180K) making it a leading quantum spin liq- 
uid candidate^ - — Vesignieite also possesses an isotropic 
kagome lattice and exhibits similar behavior to herbert- 
smithite over a wide temperature range (J ~ 53K).— ~— 
However, it was recently observed that vesignieite ex- 
hibits an exotic phase at very low temperatures which 
appears to be described by the coexistence of dynamical 
and small frozen moments^ The DM interaction, with 
an energy scale set by D, may play an important role in 
understanding the low temperature phase in vesignieite 
(D/J « 0.14). 63 Volborthitc is characterized by a small 
spatial anisotropy in the exchange interactions between 
spins, and no magnetic order has been reported down to 
temperatures of order 2K (with J ~ 77K)j£^— However, 
a signature of frozen magnetic moments appears below 
this temperature scaled - — The low-temperature mag- 
netic order may be a consequence of spatial anisotropy 
in the exchange interactions i 36 ' 70 One of the main ad- 
vantages in studying volborthite over herbertsmithite 
and vesignieite is that it can be prepared with fewer 
impurities^ - — In another deformed kagome compound, 
Rb2Cu3SnFi2r^ a valence bond state of a pinwheel type 
was experimentally observed at low temperatures 1.7 K, 
and again no magnetic order was seen over a wide range 
of higher temperatures (J ~ 180K with D/J rj 0.2)22. 
Finally, the material Cs2Cu3CeFi2 is a buckled kagome 
lattice with a quasi-one-dimensional structure in the ex- 
change constants (similar to that shown in Fig. Q} which 
shows clear signs of order at low temperature. 71 Taken 
together, these data indicate that while these systems 
may appear to be heading towards quantum spin liquid 
states over a wide temperature range, weak interactions, 
such as spin-orbit coupling, can eventually select an or- 
dered phased It is important to note that in real ma- 
terials (which are three dimensional) the DM interaction 
may take a slightly different form than the strictly two- 
dimensional form we study here. However, provided the 
layers are sufficiently weakly coupled on the energy scales 
for which we find ordering, our results should remain un- 
changed. 



B. Overview of theoretical approach and summary 
of results 

In this paper we focus on the spatially anisotropic 
spin- 1/2 kagome lattice model. Previous studies in this 
limit have neglected the DM interaction! 36 ' 73- — In our 
work, we closely follow th e q uasi one-dimensional ap- 
proach developed in Refs. |33l - l3q ] based on the appli- 
cation of the powerful one-dimensional (non-Abelian) 
bosonization techniquej 38 ' 39 In this approach one con- 
structs a two-dimensional system by weakly coupling one- 
dimensional spin chains together. Thus, the exchange 
interaction between spin chains is assumed smaller than 
the exchange interaction between spins on a given spin 




2i-l 2i 2i+l 



FIG. 1. (color online) The kagome lattice with spa- 
tially anisotropic exchange couplings. Sites along the one- 
dimensional spin chains (running horizontally) are labeled by 
i as indicated, and different horizontal chains are labeled by 
y. The scale of the antiferromagnetic nearest neighbor spin 
exchange along the chains is J, and the scale of the corre- 
sponding exchange between spins on a chain and the "mid- 
dle spins" (between chains) is J' as shown. We consider the 
anisotropic limit J' <C J in this paper. Our convention for 
the staggered Dzyaloshinskii-Moriya (DM) interaction is indi- 
cated by the different colors on the " up" and " down" triangles 
of the kagome lattice, with the up triangles having a DM vec- 
tor pointing out of the kagome plane and a down triangle 
having a DM vector pointing into the plane of the kagome. 
The arrows indicated the direction going from site i—tj. The 
strength of the DM interaction is given by lf?Jl. and shares the 
same anisotropy as J and J' (i.e. D z <C D z ) due to the spa- 
tial anisotropy of the lattice. Black arrows and brown dashed 
rectangles correspond to a DM interaction driven spin order- 
ing we found in this paper. Arrows stand for the spin direction 
while brown dashed rectangles describe dimers. Shown is the 
"Neel+dimer" order given in Eq. (|51|) . 



chain. (See FigflJ) In a one-dimensional spin-1/2 chain, 
the spin degrees of freedom can be represented in the 
spin-current formalism (non-Abelian bosonization) and 
the low-energy theory of the spin chains is described 
by a SU(2) 1 Wess-Zumino-Novikov-Witten (WZNW) 
theory^ The inter-chain interactions can then be ana- 
lyzed by perturbation theory using a scaling and a renor- 
malization group approach. The basic structure that 
emerges from this analysis is the following) 3 - 6 , there are 
two temperature scales that enter into the magnetic or- 
dering. First, the middle spins order in either a spiral 
state with small wave vector or a ferromagnetic state 
(spiral state with vanishing ordering vector) on a temper- 
ature scale T m ~ (J') 2 / J. As the temperature decreases 
further the spin chains also begin to order on a scale 
T c h ~ [J') 4 / J 3 , with an order dictated by the magnetic 
state of the middle spins^ Because of this separation 
of energy scales, the middle spins produce an effective 
exchange field on the spins of the spin chains which ulti- 
mately plays an important role in the magnetic order at 
the lowest temperatures. We note that we do not find a 
spin-liquid ground state in any regime we study. 

Our main results are the following. For weak spin- 
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orbit coupling, the DM interaction selects the plane of 
the kagome lattice in which the middle spins take on a 
spiral order (the plane of the kagome lattice itself). Our 
analysis of the spin chains shows that the spins along 
a chain are predominantly anti-parallel to the direction 
of the middle spins. However, there is a small compo- 
nent which is perpendicular to the main order, and this 
component undergoes a phase transition as the DM inter- 
action is increased. For small values of the DM interac- 
tion the perpendicular component orders antiferromag- 
netically perpendicular to the plane of the spiral, while 
for larger values this component becomes a mixture of 
coexisting antiferromagnetic and dimerized states and 
rotates into the plane of the spiral. We estimate the 
value of the critical DM interaction [see Eq. ([50| ] to be 
D/J ~ {J'/Jf/ 2 . 

The remainder of the paper is organized as follows. In 
Sec |TI] we begin with a tight-binding model with spin- 
orbit coupling on the kagome lattice. We then derive 
the effective spin Hamiltonian in the limit of a large on- 
site repulsion at half-filling. Thus, the parameters of 
a spatially anisotropic kagome antiferromagnet are de- 
termined. In this section, a non-Abelian bosonization 
mapping is introduced that we will use throughout the 
remainder of the paper. In Sec. IIIII we perform per- 
turbation theory in the exchange interaction J' between 
spin chains and middle spins. Our approach closely fol- 
lows that of Ref . [3(| . In Sec. IIVI the order of the spin 
chains and the two dimensional system is analyzed from 
the point-of-view of perturbative renormalization group 
equations. Finally, we present our main conclusions in 
Sec.[V] Some technical formulas are given in Appendix El 
and Appendix iBl 



II. MODEL HAMILTONIAN 

In order to remind the reader of the microscopic origin 
of the DM interaction, and to establish some relationship 
with other phases where spin-orbit coupling plays a role 
(such as topological band insulators^—), we first consider 
a tight-binding model on the spatially isotropic kagome 
lattice with intrinsic spin-orbit coupling. We study a 
model with only on-site interactions, and assuming half- 
filling we derive an effective spin Hamiltonian. The spin- 
orbit coupling translates into a DM interaction which 
induces an anisotropy in the exchange interaction. With 
insight from the spatially isotropic case, we then describe 
the case of a spatially anisotropic kagome lattice, which 
is the main focus of this paper. 



A. Tight-binding model on the kagome lattice with 
spin-orbit coupling 

In this section we give a derivation of the nearest- 
neighbor spin Hamiltonian in the presence of spin orbit- 
coupling on the kagome lattice. The calculation is 



standard ) 77 ' 81 but clearly illustrates how the DM inter- 
action is obtained and gives the precise form we will use 
in this work. 

For a one-electron tight-binding model on the kagome 
lattice the nearest neighbor hopping Hamiltonian is given 
by 



Hi 



etc,' 



{ij)a 



(1) 



where t is the hopping, and Ci a is an operator that an- 
nihilates a fermion on a site i with spin a. Likewise, c\ a 
is the analogous fermion creation operator. The nota- 
tion (..) denotes nearest neighbor sites, and will be used 
throughout the paper. 

The symmetry of the kagome lattice allows a nearest- 
neighbor spin-orbit hopping term of the form 



H so — i\ s 



(ij)aP 



(2) 



where X so parameterized the spin-orbit coupling 
strength, and = ± depending on the direction of 
hopping- "+" when the third site in a triangle is "to the 
left" and "-" when the third site in a triangle is "to the 
right" (see Fig. [TJ) . 78 iZ2, Here a z is the Pauli spin matrix 
representing the ^-component of the spin. 

We are interested in the strongly interacting limit at 
half-filling, which is a Mott insulator. We model the 
interactions by a Hubbard on-site repulsive potential 



H U = U'£ 



(3) 



where rii a = c\ a Ci a . In the limit of t <C U, we perform 
perturbation theory in t and obtain the spin Hamiltonian 
with anisotropic spin exchange (even in the present limit 
of an isotropic lattice) , 

Hm = JE {13) [cos(20)(Sf S? + S?S*) + SfS* 

- i / 1J sk(M)[S 1 xS J ] 2 ], (4) 

where 

4(< 2 + AL) 



J 



and 



U 



arctan(A so /t). 



(5) 



(6) 



The first line in Eq.Q is the weakly anisotropic (when 
A so <C t) exchange coupling between nearest-neighbor 
spins, and the second line is a staggered DM interaction 
with strength, 



D z = 



-Jvij sin(20). 



(7) 



In this paper we are going to focus only on this DM 
interaction, which is an intrinsic property of the two- 
dimensional kagome lattice. We would like to remark 
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that another way to obtain the results described above is 
to note that the spin-orbit coupling can be absorbed in to 
a definition of the first-nearest neighbor hopping matrix 
element t, making it spin-dependent.-£i The equations 
(gl)-© describe the spin Hamiltonian of the spatially 
isotropic kagome lattice. When the lattice is anisotropic, 
the form is slightly modified as we describe below. 



B. Spatially anisotropic case 

In the previous section we derived the Hamiltonian of 
spins on the isotropic kagome lattice, and found the DM 
conventions are of the same form earlier studied in the 
literature4^— Let us now describe how the interactions 
are modified in the spatially anisotropic case which will 
be studied in the remainder of the paper. 

In a spatially anisotropic lattice such as that shown 
in Fig. [TJ the bonds between horizontal chains are 
"stretched" leading to a modified hopping t' < t (and 
spin-orbit coupling) along the diagonals. (In fact, the 
spin-orbit coupling will also be modified along the chains 
as well because its value is dependent on the position of 
the "middle spins".) The modified hopping parameters 
and spin-orbit coupling will change the exchange con- 
stants in (|5|) and ([TJ) along the diagonal bonds. We take 
the exchange between spins on diagonals to be J' and 
on horizontal lines to be J, as before. The DM interac- 
tion also takes on a similar anisotropy with an angle 9 
described as in (slightly modified from the isotropic 
case since the spin-orbit coupling is modified, but this 
will not play a crucial role in our study) along the chains, 
and a similarly defined 9' (in terms of tf and X' so ) along 
diagonals, giving 

Dfj = —J'vij sin(20'). (8) 

Thus, in the spatially anisotropic case the Hamiltonian 
has the same form as ([3]), only with different exchange 
and DM interactions along chains and between spins on 



the chains and the middle spins. Since we are interested 
in the limit J' -C J, we will first focus on the Hamiltonian 
of the spin chains and later treat the interaction with 
middle spins perturbatively. The spin chain Hamiltonian 
is thus 

#o = JE m , v [cos(26)(S? <y Sl y + Sf >y Sl y ) + S* y Sl v ] 
+ Jsm(2e)J2 {ijhy vV.[ Sity xS j , y } z , (9) 

with J > given by ([5]) as before, the index y labels 
the different spin chains, and vf- = ±1 which alternates 
from bond to bond. For even and odd spin chains the ufj 
differs by a sign for the same (ij) bond. These features 
are evident from Fig. [TJ 

We now absorb the DM term into the first term of © 
by rotating spins with a local spin rotation^ 

/ cos(0) - sin(0) \ 
Rdm(0) = sin(0) cos(#) , (10) 

A °. v 

which rotates the spin on site i on chain y as S^j, — > 
Rdm ((-iy +v+1 9) S ilV . After the rotation, the Hamil- 
tonian © becomes an isotropic Heisenberg model 

H Q = J &i,V ' &j,y> (11) 

(ij) ,y 

which is very convenient for performing perturbation the- 
ory in the interactions with middle spins. For a purely 
one dimensional system, this result of absorbing a stag- 
gered DM interaction was obtained in Ref . [8jJ , which also 
emphasized that the DM interaction can play an essential 
role in spin ordering when frustration is present. As we 
now show, that general result also applies to the model 
we study. 

We start with the interactions between the spin chains 
with middle spins (those lying along the diagonals in 
Fig. [U between the spin chains). These terms consist 
of XXZ-type terms and the DM interaction. Performing 
the rotation (TTDl on the spin chains as described above, 
we arrive at the following form of the interaction 



V = J' COS(X) £ i>y 4^1/2 M /2 



s: 



2i,2y^l 



s: 



2i±l,2i/=Fl 



c x (y) , q x (y) 

°2i,2y ' °2i±l,2y 



i,y °2i±l/2,2y^l/2 l u 2i,2y^fl 



[si 



■ l 2i±l,2yT 



1 + Sli2y + 



s, 



2i±\,2y\ 



J' Sin(x) Yli,y [~[ s 2i±1/2,2jj=f1/2 X ^>2i.2y^l]z + [s2i±l/2,2j/Tl/2 X ^2i±l,2y^fl]i 
+ [ s 2i±l/2,2y=Fl/2 X ^2i.2y]z ~ [ s 2i±l/2,2y^fl/2 x S2i±l,2y]z] , 



(12) 



where \ = 2$' + 9 and the small "s" describes the spin degrees of freedom of the middle spins. Note that the first 
two terms describe the nearly isotropic exchange interaction (for weak spin-orbit coupling), while the last two terms 
describe the DM interaction between middle spins and the spins on the spin chains. 

Focusing on the case, J' ^ J, we can apply powerful bosonization methods to describe the low-energy physics of 
spin chains and construct perturbation theory in the interaction inter-chain interactions. The basic approach has been 
successfully applied before to anisotropic systems without spin-orbit coupling.— We begin by taking the continuum 
limit of local spin degrees of freedom which leads to a sum of a uniform and staggered magnetization^ 



a Q [M(x) + (-rfN(x)] 



(13) 
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where x = iao, with ao being the lattice spacing, and i labels sites along the chains. The definitions of these operators 
and the fusion rules which they obey are heavily discussed in the literature^ 3 . - — For completeness, we provide the 
necessary details in Appendix lAl In this representation the Hamiltonian (fTT|) is described by the WZNW SU(2)i 
theory. 39 Substituting the low-energy form (|13l) into ([12"]) we find, 



Vi = 7i E x , y 41/2,2^1/2 (m^ i{ 2x) + M$\2x)) + 7i E x , y ^±1/2,2^/2 (Mi yTl (2x) + Mi y (2x)) 
V2 = 72 £*,„ [-^ ± i/2 l2 ^i/ 2 (^(2^) - Nl yTl (2x)) + s y 2x±1/2M/2 (N% y (2x) - N^(2x)) 



(14) 



^3 - 73 Y. x . y {±)sl x % 2 . 2y ^ 1/2 d x {N x 2 ^{2x) + N^(2x)) + 73 E x , y (±> z 2x±1/2 , 2 „i /2 d x (N$ y (2x) + A% Tl (2x)), 



c (y)t 



c (y) 



where the last interaction, V3 , is important as it will gen- 
erate relevant inter-chain couplings in the renormaliza- 
tion group analysis to follow^ 

In (Ti"4"l) we have used the approximations M(2x+ao) « 
M(2x) and N(2x + a ) « N(2x) + <fd x N{2x). The cou- 
pling constants in 1(13)1 . 71 = 2a J', 71 = 2ao</' cos(x), 

2 2 
72 = 2a J'sin(x), 73 = —%-J', and 73 = -^J'cos(x) 
will flow under the renormalization group. Since \ = 
28' + 9, the strength of the spin-orbit coupling terms par- 
tially set the scale of the initial coupling constants. We 
now proceed to treat the interaction V = V\ + V 2 + V3 in 
perturbation theory. 

Note that V 2 = if the DM interaction vanishes. The 
term V 2 will turn out to give rise to dimer correlations 
for sufficiently large DM interaction, and therefore may 
potentially drive a quantum phase transition. 



III. PERTURBATION THEORY AND THE 
ORDERING OF THE MIDDLE SPINS 

The remainder of the paper is devoted to an analysis 
of the phase diagram of the anisotropic kagome antifcr- 
romagnet with spin-orbit coupling based on the pertur- 
bative interactions in (fT4")l . The basic elements of our ap- 
proach closely follow that of Ref. 36] . The new element 
in the present study is the presence of a DM interaction. 

The low-energy form of the spin chain Hamiltonian 
([lip can be written using non-Abelian bosonization (spin 
current operators)^ as a Sugawara Hamiltonian with a 
backscattering term 3 ^ 

2-ku f f 
H = -y / dx [3 R ■ J R + J L ■ J L ] + g bs / dx 3 R ■ J L , 

( 15 ) 

where u — ttJciq/2, and the second term in this expres- 
sion is the backscattering with strength g\, s < 0. 36 In this 
section the backscattering will not play any role. How- 
ever, it will be important in subsequent sections. 

The operators M and N [from Eq. flTlJl) ] that enter ex- 
pression (|14p have scaling dimension 1 and 1/2, respec- 
tively. That is, under rescaling of coordinate xt — xe~ l 
and time ti = re~', where b 



r ,df 



they change as 



M y (x,r) = b' 1 M y (x/b,T/b), 
N y (x,T) = b- 1/2 N y {x/b,T/b). 



(16) 
(17) 



There is another important operator e y (x,r), the dimer- 
ization operator, related to the continuum limit of the 
scalar product of two neighboring spins, Sj l2/ • Si+i^ — > 
(—l) x e y (x), which does not appear in the interactions 
(|14p but will be generated under the renormalization 
group flow.— It has scaling dimension 1/2, so changes 
as c y (x, t) — b~ 1 / 2 e y (x/b : r/b) under a rescaling of space 
and time. All of these operators obey fusion rules (see 
Appendix \X§ which will be used in perturbation theory 
beyond leading order. 

A symmetry analysis yields the form of all the pos- 
sible terms that are expected to be generated through 
the renormalization group (RG) transformations^ The 
most relevant ones, with lowest scaling dimension, are^ 

Hn = Yj / dxN*N* +1 + lN2 ( dxN^N*^ 
H e = ^7 e J dxe y e y+1 . 



(18) 



There are also irrelevant interactions between spin chains 
which are important because they will generate relevant 
terms of the form (flBl) 



IdNi \ dxd x N y d x N y+1 
y J 

+ X>W2 f dxd x N^d x N^l 
y J 

H M = ^2 7an / dxM z y M z y+1 
y J 

+ E7M2 / dxM^M^}. 



(19) 



With ((18)) and (1191) in hand as "targets" for impor- 
tant terms to be generated by the RG, we proceed to 
apply perturbation theory to second order in interaction 
between the spin chains and middle spins, 



Z 



-S 



1 



drV(r) + ^T j d n dT 2 V[n )\'(r 2 ) 



(20) 

where V = V\ + V 2 + V 3 is given by (JHJ). In (|2"0]1 £n is 
the fixed-point action of spin chains, and T is the time 
ordering operator. 
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At second order in the interaction, the following corn- 
s' 1 s b 

%,V 3, 3/ 



binations of middle spins will occur: Ts? A.,,. For the 



same site (i 



lzab_ 



y') the relationship Ts a {r\)s {r^) 



[0(ri-r 2 )-0(T 2 -T 1 )]e abc S c ((T 1 +r 2 )/2)is obeyed, 
where e abc is the fully antisymmetric tensor, and 6(j) is 
the heaviside function: 9(r) = 1 for r > and 9(t) = 
for r < 0. To derive the inter-chain operators one has to 
pick middle spins on the same sites with the same compo- 
nents (so that the middle spin degree of freedom "drops 
out"). Integration over the relative time is performed 
over the interval a < \tx — t%\ < ba with a — clq/u. 

To obtain contributions to the interaction V in (fT4")l 
one must select middle spins on the same site but with 
different components. After straightforward but some- 
what tedious calculations, we obtain the RG equations 



<*7i _ 7l7l _ 
dl 2iru ' c 

±13 _ _j 
dl 2 '3 n 

±h- = <y 



di 

djNl _ 
di ~ 

d*lM2 _ iMllM 2 
dl i-KU 



±12. = 1-7, 

dl 2 ' 2 



7i 

27TM ' 

±13 = -1~J„ + 

<# 2 '3 ~ 

(7awl7Ml+ 2 70W27M2) 

87ru 3 a 2 



7173 
7m 



_ 7172 
7173 

7TM ' 



(21) 



7iV2 + ii 



7JV1 

4m ' 



dl l9N2 iu , 



7M17SW2 

87m 3 a 2 ' 

7M278/V 1 

87TM 3 a 2 5 

d'YMl _ 7M2 
3? 47TM 



7j_ 

4m' 



±f8Nl 



■ ~ IS' 



which are the generalization of the corresponding results 
in Eq.(14) of Ref . (36| . (We find a small difference in the 
numerical factors of the second order terms compared to 
Ref. [HI], but they do not affect any of our overall con- 
clusions or those of that work. An example derivation is 
given in Appendix B.) 

We also obtained RG equations for exchange coupling 
constants of interaction between middle spins. We are 
not showing them here since they resemble those given 
in Ref. [3|| . The only modification due to DM interaction 
is the parametric enhancement of exchange in x — y plane. 

While the equations (|2"TT) are somewhat more involved 
than the case of vanishing DM interaction, they may be 
analyzed in a similar fashion. Because the DM inter- 
action is expected to be small compared to J', J (or at 
most of the same order), the hierarchy of energy scales 
at higher energies is not expected to change. Namely, 
there an ordering temperature for the middle spins T m oc 
(J') 2 / J, which should be relatively insensitive to the DM 
interactions. The spin chains will start to order only at a 
lower temperature T c h oc (J') 4 /J 3 ^ The procedure for 
determining the low-temperature magnetic order (that 
is, below T c h) is to first find the order of the middle 
spins. These middle spins will then produce an average 
exchange field that acts on the spin chains. Thus, the 
final ordering of the anisotropic kagome lattice will be 
obtained by analyzing the spin chains in the presence of 
a (potentially) spatially varying exchange field. 

At the energy scale for which the middle spins interac- 
tions become non-perturbative, corresponding to I ~ 1, 
one can estimate the values of the relevant coupling con- 



stants of inter-chain interactions, 

7,~"t1Sf( 1 + 2cos 4 (x)), 7*i ~ ft^cos 2 ( X ) 
7*2 ~ ftSf cos 2 ( X ) + "-^P 1 sin 2 ( X ), (22) 



■ft' Tawi 



2a. 

4u : 



where we have used that 73/vi 

- 2 2 
7mi ~ — 4 i 7M2 ~ — 4^, and substituted in the corre- 
sponding "initial" values of 71, 71, 72, 73, 73 given below 
Eq.CSD. 

In order to determine the order of the middle spins 
(which form a triangular lattice) we follow arguments of 
Ref. 36] which uses perturbation theory to derive the ef- 
fective Hamiltonian of interactions between middle spins. 
Before we proceed, however, one comment is in order. A 
full solution of the effect of the DM interactions on the 
middle spins is involved. As we have argued earlier, the 
DM interaction is not expected to have a significant ef- 
fect on the ordering temperature, T m . However, it still 
plays a role: It selects the plane of the spiral order found 
in Ref. 36] to be the x — y plane. (Without the DM inter- 
action, the plane of the spiral order is arbitrary.) With 
the plane of the spiral order determined, we then follow 
Ref. 36] and write the Hamiltonian of the middle spins 
as 



H 



A =2(J') 2 (A(1)-B(l))£ fo) 
4(J'^ 2 



x(y) x(y) 



(A(2)-B(2))Y, [if s 



x (y) a x (y) 



(23) 



where notation (..) denotes first nearest neighbor, and [..] 
second nearest neighbor pairs of spins. The coefficients 
of these terms involve 

A[r) = cos 2 ( X ) (2G M (r) + G M (r + 1) + G M (r ~ 1)) , 
B{r) = sin 2 (x) (2G M (r) - G M {r + 1) - G M {r - 1)) , 

where 



G M (r) = - 



doj 



dq Si^'f^ 1 , (24) 



and S(q,u}') is the dynamical structure factor. With the 
definition (l24l). we find 



A(r) - B(r) = f J °° dJ ft dq cos 2 (g/2) cos(gr) 

- /o~ du/ ft dq cos^r) (25) 

The first inte gral in (I25|) was numerically estimated in 
the work Ref. [361) . and it scales as 1/J. At zero tem- 
perature, the second integral in (|25j) diverges, but for 
finite temperatures small compared to J, this integral 
scales as 1/Tj 82 ' 83 Thus, overall one has that the ex- 
change coupling in ([23]), J eff = (J' 2 ) \f + C2Si ^ (x) ~|, 

where ci,C2 are order one constants. Because the mid- 
dle spins order at a finite temperature given approxi- 
mately by T m ~ (J') 2 / J, the ratio of the second in- 
tegral to the first integral in (|25p can be estimated as 
Jsin 2 (x)/T m ~ (Jsin(x)/ J') 2 , which we assume to be 
smaller than one. 
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The result of this analysis is that for temperatures of 
order T rj . it ~ Jsin 2 (x) and larger, we may use the results 
of Ref. (36J for the ordering of middle spins as a good 
approximation. Classical considerations and an analysis 
of the two-spinon dynamical structure factor suggest that 
spiral order is a leading candidate^ 

(s x ) = s [e x cos(qx) + e y sin(gx)], (26) 

where e%r y ) is a unit vector along the x(y) coordinate, 
and with q <C 1, and q — (ferromagnetic) a distinct 
possibility. Here so 1/2 is the local static moment of 
the middle spins. Below we study the effect of these two 
possible middle spin orders on the ordering of the spin 
chains: (i) spiral order in x — y plane, and (ii) ferromag- 
netic order along x direction. 

We also showed that for small temperatures, when the 
second term in J e ff dominates, the results of Ref. (3(| 
still hold. Namely, the two-spinon dynamical structure 
factor approximation still predicts the spiral order of mid- 
dle spins. 



IV. ORDERING OF THE SPIN CHAINS IN 
RESPONSE TO THE ORDERING OF THE 
MIDDLE SPINS 

We discuss how the spin chains order for q finite and 
small, and q = in the middle spin ordering described 
in Eq. (|2"6"|) . At small Dzyaloshinskii-Moriya interaction 
strength we expect a spiral order (|26l) with q <C 1, since 
the results of Ref. [36[ remain unchanged above T cr n , and 
our estimates suggest that they hold below the T crit . The 
ordered middle spins then produce an effective exchange 
field (which can be treated in mean-field theory) that 
acts on the spin chains and causes them to order. We 
now determine the order of the spin chains on various 
temperature scales as a function of the Dzyaloshinskii- 
Moriya interaction. 

To do so, we plug (|2"6"|) back into (TT2"]) and obtain the 
result that the order of the middle spins appears as an 
effective magnetic field for the spins on the chain, 

yodd/even = J> h, C0S (x) £i [cOs{qx)Sf + Sm(qx)Sf] 

± J'h S m(x) Ei(-l) 4 [cos(qx)Sf - sm(qx)Sf] , (27) 

where h = 2sq cos(g/2) and the x in the argument of the 
sin and cos is x = iao, with oq the spacing of the spins 
along the chains. The first line of (|2"T1) describes a uni- 
form spiral magnetic field while the second line describes 
a staggered, slowly rotating (in space) magnetic field. 
Performing a local rotation of spins as Sj — > R x (x)Si, 
where 

(cos(qx) sin(qx) \ 
— sin(ga;) cos(qx) , (28) 
ij 

the Hamiltonian H odd = H + V odd of the odd spin chains 



becomes, 

H odd = Jcos(q) J2 t S, • S m + J(l - cos(g)) £, S?Sf +1 
+ Jsin(g)£ i (SVSf +1 -SfSV +l ) 
+ J'hE l Mx)Sf + (-1)* sin( X )Sf) , (29) 

where the first line describes anisotropic exchange be- 
tween spins, the second line describes an effective uni- 
form Dzyaloshinskii-Moriya interaction, and the third 
line are uniform and staggered magnetic fields resulting 
from exchange interactions with the middle spins. For 
even chains, one may take x ~^ ~X- Note that since the 
effect of the spin-orbit coupling is included inx — 29 + 9' , 
the main role of the DM interaction is to give rise to the 
staggered magnetic field in the last line of (|29p . 

For the convenience of the following analysis we make 
a 7r/2 rotation about the y-axis which transforms Sf —> 
Sf and Sf — > —Sf. We further assume that q 1 
which allows us to neglect anisotropy J(l — cos(g)) up to 
corrections of order q 2 . 

We now bosonize the spin chain spins in the non- 
Abelian spin-current representation. The Hamiltonian 
(fT5f of spin chains under the rotation (|28|) has the same 
form but with updated constants: u —> cos(q)u, and 
gbs — > cos(q)gb s . We note that the backscattering cou- 
pling constant is negative, gb s < 0. The effective (uni- 
form) Dzyaloshinskii-Moriya interaction (after the rota- 
tion about the y-axis) in (USD takes the formSi 

H' = Jsin(g) J2 (SfSi+i ~ ^ £f +1 ) = dj dx(r R -J x L ), 

(30) 

where d = Jsin(g)^. The uniform magnetic field (after 
the rotation about the y-axis) in |29|) becomes 

H" = fhcos( X ) J2 S i=~ h l dx ( J R + JZ l)' ( 31 ) 

i 

where h = J'hcos(x), and the staggered magnetic field 
takes the form 

H st = J'hsm(x)J2(-iySi =lnj dxNy, (32) 

i 

where j n — J'hsin(x)- The corresponding Hamiltonian 
for the even spin chains is obtained by simply putting 
X — > — X in expressions (j3Tj) and ((321 above. 

We now perform a chiral rotation of the spin currents 
to absorb the Dzyaloshinskii-Moriya term (|30p into the 
uniform magnetic field term, pip . The right and left spin 
currents can be rotated independently while keeping the 
WZNW field theory invariant. We make the rotation 
Jfl -> RrJr and 3 L -> R l Jl using 

(cos(0) +sm(<f>) \ 
10, (33) 
±sin(0) cos(0) J 
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where <f> — arctan(d//i). Under this transformation the 
uniform Dzyaloshinksii-Moriya interaction and uniform 
magnetic field combine to take the form 



H' + H" = <Jh 2 + d 2 / dx{J z R + J 



'L)i 



(34) 



and for convenience we define h z — Vh 2 + d 2 . The 
backscattering term transforms as 

H bs = g bs J dx (i(cos(20) - 1){J%J$ + J R Jl) 
+ i(cos(20) + l)(J it J+ + J+Jl) 

+ cos(20) jyi), 
where J^ 1 = J x ±iJ y , and we define 

.91 = -(cos(20) - l)g bs , 



g 2 = cos(2(f>)g bs , 



(35) 

(36) 
(37) 



whose flow under the renormalization group will turn 
out to be important in the eventual analysis of the low- 
temperature phase of the spin chains. The fields with 
scaling dimension 1/2 transform as^£ 

N x & -+N x( - z) , (38) 
N y -4 cos(<f>)N y + sin(0)e, (39) 
e -> cos((/>)e~sm(<j))N y , (40) 

so that the staggered magnetic field becomes 

H st = 7 „ J dx (cos(cj))N y + sin(^)e) . (41) 

The inter-chain interactions transform as 

H N + H c =jdx ( lN1 N z N z +1 + j N2 N*N* +1 

+(7iV2 cos 2 (0) + 7e sin 2 (0))jV|A^ +1 

+{lN2 sin 2 (0) + 7 e cos 2 (<j)))e y € y+1 

+ cos(0)sin(,/))( 7w2 - 7e)(^2fc„+i + evJVjW) • (42) 

The total Hamiltonian, which is our subject for further 
analysis, is a sum of (|3i|). (1351) . (|4T|). (|42|) , together with 
the free part expressed through spin-currents f)15|) . In 
particular, the staggered magnetic field, Eq.(j4"Tj). plays 
a key role in our result Eq. ((5T|) . We now absorb the 
effective magnetic field (|34"|) into the spin-currents by a 



shift of the bosonic field^ <p s — > ip s H 
currents are then transformed as 

t± _ t± ^ixh'/u 
J R ~ J R C ' 

t± t± ±ixh z /u 

,J L ,J L e i 

47TW 



/2ttu ' 



The spin 



(43) 



The fields with scaling dimension 1/2 transform as2£ 

N z = cos (xh z /u) N z - sin (xh z /u) e, 
e = cos (xh z /u) e + sin (xh z /u) N z , 



and iV ;r ' J ' remain unchanged. This shift renders the first 
and third terms in ([33]) oscillatory on scales x > u/h z so 
that they will not contribute to the long distance, low- 
energy flows of the coupling constants. 



Analysis of the low-energy physics of the spin 
chains 



We will now use perturbation theory to study the cou- 
pling of the spin chains to the middle spins. We first 
note that there are two scales in the problem. At small 
length scales t < t , where t ~ ln(J/J') which is ob- 
tained as aoe e h z /u ~ 1, oscillating factors do not play 
a role. One may assume that h z = in the relations 
given immediately above. The coupling constants under 
renormalization in this case are given by 



di 



9bs 

2-7TU ' 



<ttn _ (3 

di \1 



MjL.) ry 

4ttu/ 



*7iVi(2) /i g\, s \ 

Si = I 1 ~ 4^)7iVl(2), 

*Z« _ (i 



4ttu i > ei 



and should be stopped at ^ 
backscattering we get 

9bs(l) 



.V/;.s 



1 - 



(45) 

Solving for the 
(46) 



2-jtu 



which is a decreasing function of I since g\ a < The 
solutions for the coupling constants at I ~ Z* are 



7iVl(2) = 7at 



o r / 

JVl(2) e ^ 



r _o \ 1/2 
2-ku 



* / a* o \ —3/2 
T f = I 1 - ^ 

= o ar/2 A _ rsLV /2 



(47) 



The coupling constant "f n has the largest scaling dimen- 
sion 3/2 [from (|32|) ]. We thus expect that it will grow at 
long length scales and compete with the inter-chain cou- 
pling constants, 7^1(2) (which have scaling dimension 1), 
and result in a phase transition. 

To study the phase transition as a function of the 
Dzyaloshinskii-Moriya coupling we analyze renormaliza- 
tion group flows at large scales. At large scales £ > I* 
rapid oscillations of various terms in the Hamiltonian [the 
first and third terms in (|3"5j) ] causes them to average to 
zero and we neglect them. The inter-chain interaction 
(|4"2")l after the rotation (|4~4"|) becomes 



H t = Jdx ( lx N*N* +1 + ly NIN, : 



v 

y+i 



(48) 



(44) 



where we have defined the coupling constants as "f x = 
7^2, ly = 7Ar 2 cos 2 (0)+7 e sin 2 ((/>), 7+ = 5(7^2 sin 2 (<?!>) + 
7 e cos 2 (</>) + 7jvi)- We only keep g\ and g 2 terms in the 
backscattering (13"51) . and only the first term, with N v , in 
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staggered magnetic field (|4Tj) and label its coupling con- 
stant 7„ = cos((/>)7„. The remaining terms are oscillatory 
and therefore ignored for this long-distance analysis. 

The resulting renormalization group equations describ- 
ing the coupling constants of the non-oscillatory terms 
are 



91 92 
2ttu ' 



M 



2txu ' 



(49) 



±Ul -/l± -23- - -22- ~ Hu. — (\ — _2l_ _ _22_"\ ~ 
dl \ ~ 2nu Attu) < x i dl \ x 2ttu iiru ) 'V 



d'y. 



^ = (1 



iTTU ) 



d"y T1 



'V 



-22-) -7 



The analysis^ of these equations show that there are two 
competing operators: j x and j n . The latter is due to the 
Dzyaloshinkii-Moriya interaction (recall j n — cos(</>)7„ — 
cos(</>)J'2s cos(<j/2)sin(x) and j x = 7^2 ~ (J') 4 /J 3 
are the bare couplings before the short-scale renormal- 
ization). 

From the short-scale renormalization, new "bare" cou- 
plings emerge from which the long-scale renormalization 
flows begin. These are determined from (|47| . We give a 
rough estimate of the value of Dzyaloshinskii-Moriya in- 
teraction constant at which the phase transition between 

j x and 7„ dominated states occurs. Comparing lengths 

1/2 

at which 7„e 3f / 2 ~ 1 where j n ~ sin(x) (-77) from 

(|47|). and j x e i ~ 1 where j x ~ we get the critical 

value of the angle 9 which defines the critical strength of 
Dzyaloshinkii-Moriya interaction in our model: 

/ V \ 5/2 

sin( Xc ) = - , (50) 



where for simplicity one may assume that 9^9' since 
they are parametrically of similar order of magnitude, 
giving Xc ~ 39 c . The condition (f50j) is one of the central 
results of this work. For 9 < 9 C coupling constant j x 
dominates, while for 9 > 9 C staggered magnetic field 7„ 
dominates. 

A rough numerical estimate of 9 C can be made since 
sin(26>) « D/J [from Q], which gives sin(36> c ) « 
ZD/{2J) w (J'/J) 5 / 2 . For D/J « 0.2 (possible in the 
class of materials described in the introduction), one finds 
the critical point would be reached for lattice sufficiently 
anisotropic that J'/J < 0.7. While in experiment the 
precise value of J'/J is not presently known for various 
anisotropic materials, it is possible that the application 
of pressure could further distort an anisotropic lattice to 
reach these values and possibly drive a phase transition. 
If D is smaller, the critical ratio of J'/J will in turn be 
smaller. However, we expect our main conclusions to be 
qualitatively valid even for J'/J w 0.7, in spite of our 
initial assumption J' <t; J. 



B. Resulting order of the spin chains at low 
temperatures 

For 9 > 9 C , one finds from flUJ), = 
cos(4>)(— l) v Mi, where M\ ~ ^sin(36'). Tracing back 



all the transformations performed (|43[) . (|33|) . (|28|) . and 
(fTU|). we obtain the final order of the spin chains, 

(s x x ,y) = -^lcos(0)cos( g * - (-ir+ye) 



-(-l) x+y Mi cos(0) sm(qx - (-l) x+y 6), 
(S V x ,y) = -^co S (<P) S m(qx - {-if+yff) 

+{-l) x+y M 1 cos(0) cos{qx - (-l) x+y 9), 

(S X J = 0, 
(e J/ )=-sin(</,)(-l)^M 1 , 



(51) 



= where h z = as given below (|34|l . This spin- 

orbit coupling dominated phase is characterized by a 
non-zero value of dimerization operator and staggered 
magnetization (Neel order) in the plane of the spiral (as 
opposed to Neel order in the plane perpendicular to the 
plane of the spiral as occurs for 9 < 9 c -see below). The 
presence of a non-vanishing dimer correlation in (1511) im- 
plies spin-Peierls-like physics occurs simultaneously with 
Neel order. This Neel + dimer phase is a consequence of 
the completely broken spin rotational symmetry realized 
by a uniform Dzyaloshinkii-Moriya interaction, and uni- 
form plus staggered magnetic fields. Schematically, the 
obtained phase (ISTl) is depicted in the Fig. [1] 

The same phase, namely Neel + dimer, was obtained 
in double-frequency sine-Gordon model^ and in the 
work by Garate and Affleck on spin chains with uni- 
form Dzyaloshinkii-Moriya interaction, uniform Zeeman 
magnetic field and exchange anisotropy^ both strictly 
one-dimensional systems. In particular, the Neel+dimer 
phase reported in Rcf. 86] possess Neel order which is 
perpendicular to the applied field (in our case played by 
the exchange field of the middle spins) and the DM di- 
rection. These are the same features we find in our (|51l) . 
In particular, the small q limit of our (1291) is rather simi- 
lar to Eq.(2.2) of Ref . [861] . The chief physical differences 
is that our case has staggered magnetic field and stag- 
gered DM interaction acting along the chain, while the 
Hamiltonian in Ref. [86] considers a uniform field and DM 
interaction. 

In the opposite physical limit, when the coupling con- 
stant j x dominates for 9 < 9 C , we find that the order is 
given by 



-— cos(</>) cos(qx + (-l) x+v 9), 

2lTU 



(S x ,y) 
(S x ,y) = 

(S* >y ) = -(-lf + *M 2 , 



cos(4>) am(qx + {-l) x+y t 

2lTU 



(52) 



where M2 ~ ■ I R addition to the predominant 

spiral order in the rcy-plane, this phase is character- 
ized by a Neel order along the z-direction with moment 
M 2 . The result (1521 is consistent with the prediction 
of Ref. (3(| , which does not include the Dzyaloshinskii- 
Moriya interaction (vanishes in the limit 6 — > 0). These 
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Neel+dimer 
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FIG. 2. (color online) Schematic phase diagram below the 
spin-chain ordering temperature, T c h, as a function of the 
Dzyaloshinskii-Moriya interaction parameterized by the angle 
6, as given in ([6| with Xc = 3# c - "Spiral" refers to the order of 
the middle spins which is established at a temperature scale 
of T m ~ {J'f/J > T ch , while "Neel" and "Neel+dimer" 
refers to the dominant order (on top of the spiral) along the 
spin chains below a temperature T c h ~ (J 1 ) 4 / J 3 . The "Neel" 
phase is described by (|52[) . while the "Neel + dimer" phase is 
described by (|5T|) . 



results are summarized in Fig(3] The main role of the 
Dzyaloshinskii-Moriya interaction in the Neel phase is to 
choose the direction of Neel order parameter. However, 
it also enters the spiral configuration (1521) of spin chains 
via the argument of the sin and cos. 



V. CONCLUSIONS 

In summary, we have studied the spatially anisotropic 
spin-1/2 kagome antiferromagnet with a staggered 
Dzyaloshinskii-Moriya (DM) interaction as shown in 
FigQ] As a function of increasing DM interaction, there 
is a phase transition from a spiral phase in the plane of 
the kagome lattice with an out-of-plane Neel order to a 
phase with with spiral order and an in-plane Neel order 
coexisting with dimer order. Our results are summa- 
rized in Fig. [2] A rough numerical estimate the the DM 
interaction need to enter the "Neel+dimer" phase sug- 
gest that it may be in an experimentally accessible range 
(J'/ J w 0.7), and could possibly be tuned through the 
application of pressure to a spatially anisotropic system. 
We hope this work will help inspire further experimental 
studies of spatially anisotropic kagome lattices. 

It might also be interesting to realize a phase where 
both Neel (out-of-plane) and spin-Peierls orders coexist 
for an anisotropic two-dimensional lattice. It would seem 
that necessary conditions for such phase is a fully broken 
SU(2) spin rotational symmetry (possibly obtained by 
the application of an external magnetic field) as well as 
spin anisotropy-inducing terms such as DM interaction. 



We leave that possibility as a topic for future study. 
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Appendix A: Non-Abelian bosonization and Fusion 
rules 



The operators used in the paper are defined by 

M = J^ + Jl, 
N = ^Rltr^pLp + h.c, 
e=i(RlL a -h.c), (Al) 



where R and L denote fields of right- and left-moving 
fermions in one dimension. See Ref . [39| for their defi- 
nition in terms of bosonic fields. All of the operators 
defined above are functions of position and time (x,t), 
which we have not explicitly indicated. Correlation func- 
tions of these fields are given by the following expressions 

F R = -(TR a (x,T)Ri(0,0)) = - 



Fr 



-{TL a {x,r)LU0,0)} 



2-Ku{r—ix I u+aa T ) 

1 

2iTu{r J rix /u+ot(T T ) 



where a — ao/u is short time cut-off, ao is a short dis- 
tance cut-off (distance between neighboring spins on a 
lattice), a T = sgn(r), and u = nJao/2. Below we set 
y — (x, r) for further convenience. The operators obey 
the fusion rules 



Jr. 



r/MJr/M = -F R/L (y)ie" bc J R/L (Y) + $6<*F* /L (y), 
J a R / L {Vi)N h {y2) = -^F R/L (y){e abc N c (Y) + 5 ab e(Y)], 

J a R/L {yMV2) = T^F R/L (y)N a (Y), (A2) 

from which the one-loop renormalization group flows are 
obtained. Above, y — y\ — yi and Y = (j/i + V%)j2. 



Appendix B: Derivation of renormalization group 
equations 

In this appendix we show one example of a derivation of 
the second order renormalization group (RG) equations 
of a relevant inter-chain interaction, such as those given 
in Eq. (|2T|) . To obtain the RG equations one has to plug 
(fl9l in to (|2T)|) . As an example we focus on the second 
order term 
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2^7a/i79Jvi / dx l dx 2 dndT 2 M*( yi )M* +1 ( yi )d X2 N*(y 2 )d X2 N* +1 {y 2 ), (Bl) 

where, as in the previous Appendix, the coordinates 2/1(2) stand for yi >2 = (x\,ti). Using the fusion rules of current 
operators one gets 

M*( yi )d X2 N*(y 2 ) = l -e(X,t)d X2 (F R ( yi - y 2 ) - F L ( Vl - y 2 )) = ^<X,t) L- + * + j^—p^j • (B2) 
The resulting integral is then 

7A/179JV1 / dXdt e V ( X > *) e 3/+l ( X > *) /!L ^ / dT ( (x+Lry + (x-Lt)* ) 

= lMildNi 1 ^?JdXdt e y (X,t)e y+1 (X,t), (B3) 

where X — (x\ + x 2 )/2, x — x± — x 2 , t = (ri + t 2 )/2, and r = ti — T2- The integration is over r is a < |r| < 6a. One 
has this term in the renormalization group equation for 7 e in the third line of f|2 1[) . 

I 
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